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Abstract. Fabry- Perot interferometer observations of the thermospheric O I (6300 A) emission 
have been conducted from an airglow observatory at a dark field site in the southeastern Yukon Ter- 
ritory. Canada, for the period November 1991 to April 1993. The experiment, operated in an unat- 
tended. remote fashion, has resulted in a substantial data set from which mean neutral winds have 
been determined. Dependent upon geomagnetic activity, the nocturnal location ot the site is either 
cquatorward of the auroral oval or within the oval boundaries. The data set is rich enough to permit 
hourly binning of neutral winds based upon the K p geomagnetic disturbance index as well as the 
season. For cases of low geomagnetic activity the averaged, vector horizontal neutral wind exhibits 
the characteristics of a midlatitude site displaying antisunward pressure-gradient-driven winds. As 
the geomagnetic activity rises the late afternoon and evening winds slowly rotate sunward in an an- 
ticlockwise direction, initially remaining near 1(X) m/s in speed but eventually increasing to 300 
m/s for K p > 5. For the higher levels of activity the observed neutral wind How pattern resembles 
a higher-latitude polar cap pattern characterized by ion drag forcing of thermospheric neutral gases. 
In addition, rotational Coriolis loreing on the dusk side enhances the ion drag lorcing, resulting in 
dusk winds which trace out the clockwise dusk cell plasma How. On the dawn side the neutral 
winds also rotate in an anticlockwise direction as the strength of geomagnetic disturbances in- 
crease. Since the site is located at a transition latitude between the midlatitudc and the polar cap, 
the data set provides a sensitive test for general circulation models which attempt to parameterize 
the contribution of magnctospheric processes. A comparison with the vector spherical harmonic 
(VSH) model indicates several regions of poor correspondence for December solstice conditions but 
reasonable agreement for the vernal equinox. 


Introduction 

The upper thermosphere is the region which separates the 
interplanetary environment from the terrestrial domain and is 
consequently affected from above by solar-driven processes. 
Though relatively vacuous compared to the troposphere, the 
thermosphere is dense enough to support its own climate and 
weather patterns. The dynamical influence of solar heating 
tends to drive a pressure gradient wind from the subsolar point 
across the terminator towards the night side ot the Earth. At 
high latitudes, terrestrial magnetospheric processes are impor- 
tant in determining the thermospheric wind flow primarily 
through momentum transfer from rapidly moving plasma. In 
addition, the inherent Coriolis forcing associated with the ro- 
tation of the Earth has an effect. 

Long-term observational programs monitoring the neutral 
dynamics in the upper thermosphere at auroral latitudes have 
been performed directly via high-resolution optical methods or 
indirectly via incoherent scatter radar reflections. Ground- 
based optical experiments with a Fabry-Perot interferometer 
were performed at College, Alaska (invariant latitude A = 65) 
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[Sica, 1984; Sica et ai , 1986a] from which average neutral 
winds for low (A p < 1 1), moderate, and high (A }) > 39) geomag- 
netic conditions were determined. A corresponding studv hs 
Aruliah et ai [1991a] discussed average neutral winds tor dit 
ferent ^conditions and low (F {0 j ^ 100), moderate, and high 
(f, 07 > 170) solar flux levels as observed at Kiruna (A = 65 » 
Seasonal effects on the magnitude of the upper thermospheric 
neutral winds have also been investigated at Kiruna [Aruiuih n 
ai, 1991b]. In this latter study the influence of low i7 i(l • * 
110) and high (F 107 > 110) solar flux levels during different 
seasons were also studied. Mean meridional winds have been 
derived from European incoherent scatter (EISCAT) (A ~ 6"i 
radar data as a function of season during solar minimum condt 
lions 1 Titheridge, 1991 1. This set is the only study which pm 
vides summer season wind measurements. Neutral winds have 
also been measured by a Michclson interferometer at Sn 
dankyla and Troms0 (EISCAT radar sites) during short cam 
paign runs between November 1988 and January 1991 prmid 
ing similar conclusions as the Kiruna studies [Faultoi et ai . 
1993). Finally, Crickmore (1994) has presented a study ol the 
mean upper thermospheric winds at Halley (A = 62°S) as a (unc- 
tion of season, high (ZK p > 24) and low OLK p < 16) geomag 
nctic activity, and high (F U)1 > 190) and low (F107 < I 44 i vo- 
lar flux level. 

Model descriptions of the upper thermospheric neutral 
winds have appeared lor both empirical and theoretical mu. L ev 
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Hedui et al. | 1 99 \ | have collated ;i large mimbei ot satellite 
home ami ground-based observations u> generate the empmeal 
hon/onlal wind model <IIWM) description ot global iheimo 
spheric dynamics. Theoretical simulations ol thermosphei ic 
dynamics have been performed independently based upon two 
constantly improving models: the National ( enter lor Atmo- 
spheric Research (NCAR) thermospheric general circulation 
model (TGCM) | Dickinson el at.. I9SI| and the l inger its 
College London (UCL) TGCM \h tiller Howell on </ A *ee\. I9X0| 
These models are typically difficult to use. which has led in the 
development ol the VSH parametert/ed version ol the NCAR 
model {Killeen et ai, I *>87|. 

We report in this paper observations ot upper thermo 
spheric horizontal neutral winds based on a new detection 
technique, a Fabry-Perol interferometer utilizing a bare CCD 
as the detector. The technique permitted instantaneous viewing 
of multiple copies of an O I (6300 A) emission line protile 
thereby inherently increasing the signal-to -noi.se ratio over 
traditional single fringe measurements. Consequentl> . even 
though the experiment operated for only two dark observing 
seasons, a large amount of high-qualily data was collected 
These data sets provide an element towards the high latitude 
component of the coupling energetics and dynamics of atmo- 
spheric regions (CEDAR)-sponsorcd coordinated analysis ot 
the thermosphere (CAT) program [ Hagan et at., 1991 |. The re- 
sults are comparable with the auroral /one observational stud- 
ies described above, though with the greater density of data in 
the current study, it is now possible to simultaneously com- 
pare the thermospheric horizontal neutral wind as a function of 


both Aj, and season, the Inrmer discriminant being imu h finer 
than ever previously coneei vahle. 


An airglow ohservatory. shown in Figure I. was established 
m 1 99 1 at a dark held site near Watson Lake, Yukon Territory 
< latitude: 60 N. longitude: I 29 W ) . 'The instrument comple- 
ment included a Fabry -Perot interferometer, a Michelson intcr- 
Icronicter. and an , ill-sky imager, (he latter operated bv Utah 
State University. In the present study, only data acquired by 
the I ahry-Perot interferometer shall be discussed. This instru- 
ment operated in an unmanned, remotely controlled mode dur- 
ing two dark observing seasons: (I) November 1091 to April 
1992 and (2) November 1992 to April 1993. The Fabry Perot 
observation program concentrated solely on acquiring mea- 
surements of the upper thermospheric O I (6300 A) emission 
Icaturc. Much of the construction of the Fabry-Perot interfer- 
ometer experiment resembled that described by Meriwether et 
al |I9K3|, but the detection system utilized a bare charge-cou- 
pled device and a microprocessor-controlled data acquisition 
system J Niciejew ski et al.. 1994], Operational parameters of 
the Fabry-Perot interferometer arc listed in Table I. 

The CCD chosen for the detector was the Kodak K A FT 400 
chip which was supplied hermetically sealed in a thermoelcc- 
iricaliy cooled liquid circulation unit. The system ran in the 
multipinned phase (MPP) mode which reduced dark noise levels 
substantially. In addition, the individual detection elements 
i pixels) were co-added on chip to form 8-x-S-pixel-wide "super 


Kxpcrimcnt 



Figure I. Photograph ot the atiglnw ol>sci\an>i\ at the W ilson Lake field site Die F.ibiv Perot mierterom- 
elet resides on an optical bench Ivnc.iih one ot iwo ti aiispaieni I lesicl.iss donuw j hr site is located on a 
mountain lop about 10 miles ill’ bin r.i-a ol lown 
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Table I. Operational Parameters lot the Pabry-Pcrol 
Interferometer at Watson Lake Airflow Observatory 


Optics 


Parameter 

Etalon plates; reflectivity at 

6300 A* 

01 % 

Clear aperture 


12.5 cm 

Spacer gap (solid annular quail? nng) 

1 0 cm 

Typical working finesse for innermost channel 

14 

Interference filter 



center wavelength 


6300 A 

band pass 


8 A 

transmission 


70% 

Lens Diameter, cm 

Focal Length, cm 

Objective 

15 

155 

Collimating 

5 

60 

Condensing 

5 

10 


The lens combination given above implies nearly three complete free 
spectral ranges at A.6328 A imaged onto the detector. The tree spectral 
range is 0.20 A at X6300 A. The detector is a Kodak KAF1400 with 
1317 x 1035 pixels (pixel size is 6.8 x 6 8 pm): readout noise, 11.2 e~ 
RMS; dark noise, 0.003 e'/s at -42 °C 

* Manufactured by 1C Optical Systems Ltd 

pixels" prior to image readout in a successful effort to reduce 
readout noise. Using the lenses and t i Iter described in Table 1, 
three O I {6300 A) profiles, or fringes, were simultaneously 
imaged onto the detector. The integration period per sky image 
varied between 10 and 180 s based solely upon the O 1 (6300 
A) signal strength. The integration period was determined dy- 
namically on site by the data acquisition controller based on 
an initial 10-s exposure. The criterion employed was to choose 
the shortest integration which would provide a wind error of 
approximately 5 m/s but not to observe more than 180 s in 
any direction. Instrumental drift of the Fabry-Perot interferom- 
eter was monitored by frequently observing a Ne calibration 
lamp. Prior to the start of each of the two dark observing sea- 
sons, an instrumental calibration was performed with a fre- 
quency stabilized HeNe laser. Analysis ot the drift measure- 
ments and the instrumental calibrations indicated that the local 
experimental environment remained stable during each season. 

The routine patrol mode of operation of this experiment was 
to observe both the horizontal and vertical thermospheric 
wind by cycling through a series of five sky scans: a zenith-di- 
rected measurement to define the vertical wind, north- and 
south-directed measurements to define the meridional wind, and 
east- and west-directed measurements to define the zonal wind. 
The cardinal point scans were performed at a 45° elevation an- 
gle.. Following each cycle, the raw images were reduced into 
wavelength space using the technique described by Mulligan 
[1986]. This provided fringe profiles from annular or ring-sum 
bins. Geophysical parameters were extracted from the multiple 
fringe profiles using a technique based on that of Killeen and 
Hays [1984]. These parameters were then both stored in the 
data logger as well as broadcast over a telephone modem line 
for instant viewing. Cloud cover information was obtained 
from the Watson Lake airport meteorological service at 1-hour 
intervals and later used to remove poor-quality data. 

Results 

Guided in part by the CAT analysis procedures, certain bin 
sizes and conditions have been chosen lor studying the mean 


luniiral horizontal wind. To describe seasonal influences on the 
neutral winds, the "December solstice’ season was defined to 
include the 2 months December and January SimilaiK. the 
"March equinox" season was defined as data accumulated during 
March and April. Unfortunately, no data were acquired during 
the autumnal equinox period and sky conditions were too 
bright for any June solstice measurements. The Labr\ -Perot 
data sets thus provided four seasons for detailed study. Decem- 
ber solstice. 1991; March equinox, 1992; December solstice, 
1992. and March equinox, 1993. Since the solar cycle is much 
longer than the observation windows, no further filtering 
based upon the F i0 j solar flux measurements was performed. 
Table 2 defines the solar flux values for the four seasons of 
study. Since Watson Lake can be both a subauroral and an auro- 
ral site, filtering based upon the geomagnetic disturbance in- 
dex, K p , was performed. As shall be demonstrated later, the ef- 
fect of geomagnetic activity upon the neutral winds is quite no- 
ticeable. Finally, a bin width of l hour was chosen for the data 
study, based partly upon the frequency of cloud cover informa- 
tion and partly on the need for desirable statistics for the bins 

Auroral Contamination 

A potential source of error in describing the mean upper 
thermospheric winds is the contribution of auroral emission ro 
the measurement. This is a fault that is inherent to all group.,: 
based optical diagnostic techniques which attempi to obvc:\ e 
terrestrial airglow emission and is evidenced either through a 
particularly strong enhancement in the observed signal cb.m 
nel or through unwanted contributions from contaminant m 
roral-only channels, the latter summarized by Gatiinfe' 
Vallance Jones [1974], There have been several studies 
scribing an altitude change in the peak height of O I \ 

emission during aurora [5(ca el ai, 1986b, McCormac m .. 
1987a], indicating that a ground-based observation rn.m u . 
ally be viewing emission weighted to a much lower a!) ’ 
than during a "pure nightglow" case and subsequentlv p«'vr ,• * 
contaminant for the present statistical study. Consequent t 
simple technique was developed to exclude obvious jumr.i. : > ■ 
from desired nightglow-only data. 

The occurrence frequency of aurora have been stud- : ’ 
many years, and their geomagnetic time-latitude disir f 
for various levels of magnetic disturbance have been 
rized by Feldstein and Siarkov [1967]. Figure 2 -Imp 1 • 
mathematical representation of the oval for dilterem ic <■ 
the magnetic disturbance Q index (similar to the K index 
roral region stations) based on the formulae ol Hol :«> -• 
Meng [19751. Given the invariant latitude of Watson 1 .in- ‘ 
65) and that local magnetic midnight occurs at 9 65 1 1 
clear that the site is a "midlatitude" station for Q leseN ' 1 o 
and perhaps 2. At greater disturbance levels the sue hev -*> * 
the boundaries of the auroral oval. 


Table 2. Mean Solar Flux Values for the Four Seasons 
Study 


Season 

Solar Dux 


December solstice, 1991 

21 3 x 


March equinox, 1992 

164 7 


December solstice. 1992 

125 9 


March equinox, 1993 

125 H 
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ll is interesting to display the mean hourly binned () I 
(6300 A) signal levels tillered aeeoidmg to the K p index For 
this analysis and for related analyses in the remainder of this 
report, measurements were soiled into bins 1 K,, unit wide, that 
is, 0 < K p < 1. 1 < K p < 2, etc., primarily because the auroral 
oval moves equatorward over Watson Lake rapidly as geomag- 
netic activity increases. Figure 3 displays the average signal 
levels with their corresponding standard deviation values tor 
the December 1992 solstice period. There are several common 
features evident in this series. The evening twilight decay of 
the O I (6300 A) signal reflects the approximately 1 -2-hour 
time constant for the decay of the plasma densities in the F re- 
gion since, in the absence of auroral precipitation, the excited 
0(1D) state is primarily populated via dissociative recombina- 
tion of 0 2 + and e~. As a corollary, the morning twilight en- 
hancement is a much more rapid event A magnetic midnight 
enhancement in signal level is evident tor all levels ot geo- 
magnetic activity. For K p index levels greater than 2, a marked 
asymmetry in signal level about magnetic midnight exists. 
This latter behavior is likely resultant from the skewing of the 
auroral oval to the postmagnetic midnight hours as displayed 
in Figure 2. Finally, it is clear that the mean signal levels rise 
in intensity as geomagnetic activity rises, though only by fac- 
tors of about 5 rather than by orders ot magnitude. 

The behavior of the average signal level is similar for the 
higher solar flux December 1991 solstice period. The primary 
difference is that the signal level is greater in magnitude, as 
would be expected, since the plasma density which drives the 
dissociative recombination reaction is greater. An interesting 
anomaly is that the signal levels for the magnetic midnight 
periods for K p > 3 actually are decreased with respect to pre- 
and postmidnight periods. A less obvious behavior is evident 
for the corresponding period in Figure 3. In both cases it is 
possible that the equatorward edge of the oval has descended so 
far south that Watson Lake is entering the polar cap. The aver- 
age signal plots for the vernal equinox seasons are similar 
though shorter in duration. The signal enhancement at mag- 
netic midnight for \ow-K p conditions is present, as well as the 
asymmetry in signal levels about magnetic midnight. Signal 
levels also increase as geomagnetic activity increases. 

To estimate the contribution to the 6300-A signal from au- 
roral contributions, an auroral ionization (GLOW) [Solomon, 
1993] model was run. This model calculates ionization and ex- 
citation rates, energetic electron production and transport, ex- 
cited species densities, and airglow emission rates based upon 
solar and geomagnetic input parameters. The neutral atmo- 
sphere is modeled by the mass spectrometer/incoherent scatter 
(MSIS) model and the ionosphere description is obtained from 
the international reference ionosphere (IRI) model. To simu- 
late the background O I (6300 A) airglow, a run was made for 
December solstice, A p = 5 (approximately K p = 1.5) and is dis- 
played in Figure 4a. Here, the general features witnessed in the 
Watson Lake observations are present: the slow dusk intensity 
decay, the rapid dawn brightening, the increase in signal 
strength as solar flux increases, and the magnetic midnight 
enhancement observed for the low-/C^ cases. The midnight fea- 
ture arises from an increase in the IRl-modeled F region plasma 
density near the peak altitude of O I (6300 A) emission. This 


occurs ,u approximately the time that the Harang discontinuity 
passes over Watson Lake, that is, the boundary between the 
dusk and dawn plasma convection cells where the plasma flow 
is primarily and strongly antisun ward. Runs were then made for 
various simulated aurora: Figure 4b displays the resultant O 1 
(6300 A) emission for an aurora characterized with a mean en- 
ergy of 1 keV and an energy flux of 1 erg cm 2 s* 1 . As expected, 
the resultant emission is enhanced due to the auroral precipita- 
tion. in this case by more than an order of magnitude during 
dark night conditions. 

Subsequently, Figure 5 displays the mean signal levels of 
Figure 3 but in a histogram representation. Here, each hourly 
bin of reported signal measurements have been additionally 
binned according to the relative value: 0 to 5 units. 5 to 10 
units, etc., and then displayed centered upon the hour of the 
bin. The vertical scale has been normalized to the maximum 
histogram depth. As the geomagnetic index increases, the 
population of the bins slowly changes towards a bimodal dis- 
tribution with the largest proportion of contributors in the 
first two bins (0 to 5 and 5 to 10) and then an occasional minor 
peak corresponding to the higher intensity levels. For the 
purpose of excluding purported auroral O l (6300 A) contribu- 
tors to neutral wind bins, measurements corresponding to sig- 
nal levels greater than 10, a conservative value suggestive of 
the maximum airglow-dominated level, have been excluded 
from the reported averages. 

Neutral Wind Observations 

The meridional winds for the December 1992 solstice period 
are displayed in the series of plots in Figure 6. The convention 
for positive winds chosen for the display is north and east pos- 
itive. The mean winds with associated one-sigma standard de- 
viation vertical bars are displayed for each hour of available 
measurement. For geomagnetic indices greater than 4, there 
were insufficient data points during certain hourly bins; these 
are consequently left blank. In addition, as the disturbance 
level increases, the hour of peak equatorward wind shifts from 
local midnight to predawn, though the magnitude does not in- 
crease appreciably until the highest disturbance levels The 
behavior of the mean zonal winds is more complicated. As the 
disturbance level increases, the zonal wind changes from a pre- 
dominantly eastward wind to one with a more diurnal nature: 
westward in the evening and eastward in the morning. 

An alternate form for displaying the horizontal winds is to 
construct a vector wind and display such on a polar dial coordi- 
nate system. Figure 7 displays the neutral horizontal wind us- 
ing a geomagnetic latitude-magnetic local time format. In this 
format the Sun is towards the top of the panels, magnetic mid- 
night occurring at 9.65 UT at Watson Lake's position 
(invariant latitude is 65°N). The central cross is the magnetic 
pole with concentric circles surrounding it corresponding to 
latitudes 80°, 70°, and 60° respectively as the radii increase 
The wind magnitudes are normalized to the barbed arrow at the 
bottom right. Each measurement is shown with the tail located 
at 65°N and the barb pointing in the prevailing wind direction, 
that is, the azimuth angle of the horizontal wind vector. As the 
disturbance level increases, the wind direction changes suh- 


Figure 2. Fcldstein ovals for various levels of the Q magnetic disturbance index. The polar plots have been 
constructed for geomagnetic latitudes 60° to 90° and as a function of magnetic local time. Watson Lake has an 
invariant latitude of 65°, and magnetic midnight occurs at 9.65 UT. 
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Figure 3. Mean O I (6300 A) signal levels during the December 1992 solstice 
deviation is also shown with each hourly mean. Binning was performed for unitari 


period. One-sigma standard 
ly increasing K p indices. 


stantially. Initially, the evening winds are directed in an an- 
tisunward direction roughly towards 0200 MLT As geomag- 
netic activity increases, the winds slowly rotate in a counter- 
clockwise (from above) direction until, at the highest level, 
they are roughly tangential to the 65° latitude arc During the 
morning hours, the wind vectors also rotate counter-clockwise 
with their magnitude increasing substantially 

Polar dial wind constructions for the March 1992 equinox 
appear in Figure 8. The evening and premidmght winds also 


follow a counter-clockwise rotation as the magnetic distur 
bance level increases and do so to a greater extent than during 
corresponding hours of the December solstice. During the 
morning hours, the rotation is also coupled with a substantial 
increase in the magnitude of the horizontal wind. 

The neutral winds during the December 1991 period are sim- 
ilar to those shown in Figure 7. However, the general appear- 
ance of the polar dial displayed winds in Figure 9 show that the 
behavior is different by roughly I K p index unit. Note that the 
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GLOW 01(6300) signal 


Time (LT - UT 


Watson Lake 
Deo solstice 

P10 7 variation: 75. 250. 25 
Ap ■= 05 


Time (UT) 


GLOW 01(6300) signal 


Time (LT = UT - 8) 


Watson Lake 
Dec solstice 

F10.7 variation: 75. 250. 25 
Ap - 32 


Time (UT) 


Figure 4. GLOW model output of the O I (6300 A) intensity for December solstice, A (} = 5, for various solar 
flux levels ranging from F i0 7 = 75 (heavy line) to F , 0 7 = 250 in increments of 25 F ]07 units, (a) Photochem- 
istry only, no aurora; (b) with an aurora of mean energy 1 keV and energy flux of 1 erg enr 2 s 1 . 
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Figure 5. Signal level data of Figure 3 but displayed in a histogram form. Each hourly bin of Figure 3 was 
turther subdivided based on the strength of the signal level: 0 to 5. 5 to 10, etc. The resulting distributions are 
displayed normalized to the peak bin population for each K p level and ordered about the hour of the measure- 
ment. 


solar Mux values for the two solstice periods correspond to a 
solar maximum level tor late 1991 and low to moderate for late 
1992. Finally, the March 1993 horizontal winds shown in 
Figure 10 are very similar in appearance to those acquired a 
year earlier. The primary difference is that the magnitude of 
some of the wind vectors is diminished. 


Neutral Wind Model 

The vector spherical harmonic (VSH) model, mentioned ear- 
lier. was used to assess the current theoretical understanding of 
the horizontal neutral winds at a subauroral/auroral site in 
comparison to the observations. The version employed for the 
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Average meridional winds^ Average meridional winds^i^.vi 

® D*cfl 2 Oc« D«fl2 «o<«Wr 



Average meridional windsw.««r u**.n Average meridional winds Watson LoU. yT 

® 0*&2 taMjc* 0«<«2 sdstcc 




Figure 6. Neutral wind data for the December 1992 solstice period, (a) Meridional winds ordered with north 
positive; (b) zonal winds ordered with east positive. Mean wind values are shown with plus or minus one- 
sigma standard deviation vertical bars. 


comparison smoothly summarizes the output from 18 different 
NCAR TGCM computer runs by a series of spectral coefficients 
for a total of 13 different geophysical parameters as a function 
of time, position, and disturbance and activity levels. The ba- 
sis set of TGCM runs were acquired using a global grid ordered 
at 5°x5° in latitude and longitude, at 0.5 scale height in alti- 
tude, and at 5 min in time and were run until diurnally repro- 
ducible results were apparent. Runs were made for both the June 
and December solstices, as well as for a general equinox date. 
Different runs were evaluated for geomagnetic disturbance 


levels using A p indices of 5, 11, and 32 (K p values of roughls 
1.5, 2.5, and 4.25) as well as for high (F \qj = 215) and low 
(Fj o 7 = 72) solar activity cases. 

Figure 1 1 displays the VSH model neutral wind simulation 
using the geomagnetic latitude-magnetic local time coordinate 
system for December solstice conditions, F ]0 j = 125. and a 
range of A p indices. The solar flux value and the season de- 
scribe the December 1992 solstice period. The range »»i V 
indices includes the three for which TGCM runs were performed 
as well as an extrapolation to A p = 0. Though the VSH runs u* 
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Figure 6. 

quire geographic latitude-longitude inputs, the corresponding 
geomagnetic latitudes chosen lor the simulation were 55°N, 
N. 75°N, ail with a corresponding magnetic local mid- 
night equal to that at Watson Lake. The model predicts that as 
the geomagnetic disturbance level increases, the horizontal 
U)lkl vectors rotate in a counter-clockwise fashion in the noon 
■o dusk and dusk to midnight sectors. Between midnight and 
dawn the wind primarily increases in magnitude as the index 
h'cl uses, with some anticlockwise rotation nearest to 06(H) 
Ml I loi lowest latitudes. The dawn to noon sectoi (an area 
uiinh\i-i \ able from Watson Lake) indicates strengthening m 


(continued) 

the am isun ward How across the pole as activity increases, 
along with a southward extension of the polar cap influence. 

f igure 12 contains the results for the simulation corre- 
sponding to the March 1092 equinox period generated with a 
solar flux level ot 175. The wind vectors are. in general, simi- 
I a i i ° corresponding vectors in Figure 11 Concentrating 
soiciv upon the 65° circle of wind vectors, the primary differ- 
cnccN between Figures I I and 12 occur m the postdusk period 
Um the lowest A f , level and in the postmidmght hours for the 
lughci disturbance levels In parlicutai. the horizontal wind 
magnitudes between IKKX> and (MOO MU' are significantly en- 












Niai Jl WSkl hTAL.: THERMOSPHERIC HORIZONTAL NI-ITOAL WINDS 


25.S 


VSH run, F1Q.7 = 125 — Dec solstice; Ap = 0 


12 



VSH run, F10.7 — 125 — — Dec solstice: Ap - 5 

^ s 

12 



VSH run, F10.7 = L25 Dec solstice; Ap = 11 




Figure 11. Neutral wind simulation for conditions reproducing the December 1992 solstice shown in a ge- 
omagnetic latitude-magnetic local time coordinate system. The coordinates for the simulated sites corre- 
sponded to magnetic latitudes 55°N, 65°N, and 75°N, all with a local magnetic midnight similar to Watson 
Lake's. 


hanced in the equinox simulation. Finally, the behavior of the 
neutral winds at 65°N for simulations corresponding to the De- 
cember 1991 solstice {F \qj = 225) and the March 1993 
equinox (F 107 = 125) are similar to those developed for the 
December 1992 solstice and the March 1992 equinox, respec- 
tively, and are not shown in this report. As expected, the ear- 
lier solstice winds, simulated for a higher solar flux level, arc 
typically greater in magnitude, but roughly similar in direc- 
tion. Similarly, the later March equinox period corresponds to 
a lower solar flux level and consequently lower magnitude hor- 
izontal neutral winds in the simulation For December solstice 
conditions, the greatest neutral horizontal winds simulated for 


65°N occur postdusk at the highest geomagnetic disturbance 
value (A p ~ 32), are greater than 300 m/s in magnitude, and are 
directed sunward. For March equinox conditions, the period 
following local magnetic midnight contains the greatest 
winds, also of magnitude greater than 300 m/s but directed un- 
tisunward. 

Discussion 

Watson Lake resides at a transitional boundary of thermo 
spheric flow patterns, and consequently, long-term observa- 
tions of mean thermodynamic characteristics arc essential in 
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Figure 12. Neutral wind simulation for conditions reproducing the March 1 l )02 equinox period shown in a 
geomagnetic latitude-magnetic local time coordinate system. 


refining current understanding of global circulatory climatic 
models. The plasma density in the upper thermosphere is sutti- 
cient to forcefully drive neutral How patterns if solar and mag- 
netospheric processes perturb and accelerate the plasma drilts. 
The polar regions are open to the greatest plasma perturba- 
tions and through momentum transfer between collisions ot 
neutral species and ions, the polar regions are also susceptible 
to ’’weather modifications" in the neutral thermosphere. 

The current study has shown that the neutral winds at 65 'N 
magnetic latitude are highly susceptible to modification as the 
planetary geomagnetic disturbance level varies The K jt depen- 
dence in the meridional winds, shown in Figure 6. indicates 
that the hour of peak equatorward wind changes by 5 hours 
from K ft - 0 to K p = 5. The magnitude change p* not as dra- 
matic. with the peak equatorward wind remaining at 100 m/s 


between K p = 0 and K ft - 2. rising to 150 m/s from K p = 2 to K p 
= 4. and then peaking at 250 m/s between K ff = 4 and K n - 6 
As the disturbance level rises, the postdusk poleward wind in- 
creases in magnitude from 50 to 250 m/s. The zonal wind be- 
havior is much more complicated, changing from an appar- 
ently diurnal sanation at the lowest disturbance level, through 
an oscillatory behavior at middtsturbance levels, back to an 
apparently diurnal behavior at the highest K p levels though in 
antiphase with the lowest disturbance level. Consequently, it 
is not appropriate to characterize the wind components at this 
site simply by their maximum \ allies as may be done for the 
highest polar latitudinal studies \\ft Comuic et at.. I ( *87b|, 
since it is clear that the peak zonal wind magnitude actually de- 
creases with increasing distm bancc level before it increases at 
the higher K f , levels. 
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I lie VSH wind field pattern di>pla\ed tn Figuic I I fur De- 
cember solstice, a low solar to mntsulat actmiv level, illus- 
trates that two separate wind patterns are mixed at high lati- 
tudes. A predominant antisunward thus seen at equaioruard lati- 
tudes for the lowest disturbance levels is driven hv the pres- 
sure-gradient force directed from the warm sunward side to the 
colder night side of the Earth. At the highest latitudes the wind 
pattern consists of a two-eell structure split approximately 
along the noon-midnight meridian As geomagnetic activity 
increases, the flow at 65 °N changes from the pressure-gradient 
form at A p - 0 to the two-cell flow at A p -II. resembling the 
flow at 75°N for A p = 0. At this highei disturbance level, the 
flow at 75°N has changed significant!) from that at A p = 0. 
This new How at 75°N also resembles the flow at 65'’N when A p 
is increased to its highest level. To summari/e. the dominant 
flow at 65°N simulated by the VSH model changes character 
from pressure-gradient driven to ion diag driven as A p in- 
creases from 0 to 32, though the former does not truly disap- 
pear at the higher activity level. 

The degree of correspondence between the measurements in 
Figure 7 and the VSH simulation in Figure 1 1 for 65°N depends 
on the time period examined. Between 1200 and 1800 MLT. 
the observed wind direction for K p =• 0 is similar to the A p = 0 
simulation, but the observed magnitudes are significantly 
greater. In fact, the observed winds do not begin to change di- 
rection until the K p level exceeds 3. unlike the simulated 
change at A p = 5 (K p - 1.5). At the highest K p levels the ob- 
served winds do follow the simulation for A p = 32 quite well 
Between 1800 and 0000 MLT there is no correspondence be- 
tween observed winds at K p - 0 and the simulated A p = 0 winds 
at 65°N. However, the VSH simulation for 55°N strongly re- 
sembles the observed winds. The observed winds for K p levels 
greater than 1 and less than 4 somewhat resemble the simulated 
curl in the VSH winds for 65°N. Once again, the observed 
winds for the highest disturbance levels correspond well with 
the A p = 32 simulation. For the period 13000 and 0600 MLT. 
the correspondence between observations and simulations is 
quite good. In particular, the observed decrease in the wind in- 
tensity and the convergence at 0300 MLT for K p = 0 is well 
simulated in the A p = 0 case. The observed counter-clockwise 
rotation and enhancement in magnitude as K p increases is well 
matched by the VSH model runs. Apparently, the VSH simula- 
tion is best during predawn but significantly underestimates 
the inertia of the pressure-gradient force between dusk and 
midnight unless one shifts the simulation latitude equatorward. 

Regarding the March 1992 time frame, the observations in 
Figure 8 cover about 4 to 5 hours less ot the diurnal cycle than 
those during the winter solstice. However, it is sti II valuable 
to compare with the VSH output of Figure 12 excluding ihe pe- 
riod 0500 to 1600 MLT. In fact, the comparison between ob- 
servation and simulation is quite favorable to the simulation. 
The behavior in the observed winds between 1800 and 0000 
MLT for K p = 0 is more similar to the VSII run at A p = 0 lor 
55°N in terms of magnitude, but closer in direction to the 65° N 
run. As the K p level increases, the observed rotation in the 
wind vectors is predicted quite well by the equivalent VSH runs. 
For the period 0000 to 0600 MLT. the winds in, h, l . lowest dis _ 
turbance observation arc not well matched ur , 0 | q lc VSH 
runs, but as the K p level increases, the predicted winds for the 
higher A p levels match reasonably well m Jnechnn and mag- 
nitude. 

The observations in Figure 9 correspond io a high level of 
solar activity. The corresponding VSI! iun mot shown here) 


strongly resembles that of Figure II except for an increase in 
magnitude in most wind vectors. However, most ol the com 
nients i elating to the comparison of December 1003 solstice 
observations and the VSH runs for F, 0? = 125 also applv here 
In particular, between 1200 and 1800 MLT. the observed wands 
for K p = 0 arc still greater in magnitude than those simulated 
even though the directions are approximately correct. Obser- 
vations as the K p level increases do not match the 65' N VSH 
circle for higher A p levels until the highest K p levels are 
reached. The correspondence from 1800 to (XXX) MLT is much 
better for the December 1991 solstice as is the agreement with 
winds observed from 0000 to 0600 MLT. Finally, the observa- 
tions in Figure 10 match quite well with a VSH simulation tor 
F l0 i = 125 The only visible differences between the modelled 
and the observed winds are that occasionally the modelled 
winds are loo high in magnitude. 

The subtle differences for a constant geomagnetic distur- 
bance level witnessed in the horizontal wands displayed in 
Figures 7 to 10 between different seasonal representations and 
the same season but a different solar flux level indicate that the 
interplay between different forcing mechanisms is complex. In 
general, the observed trend that the horizontal winds at 65 : N 
evolve from a midlatitudinal to a polar cap flow as magnetic 
disturbances increase is repeatable for both the equinoctial and 
solstice seasons, but there are subtleties that cross the binning 
strategy employed in this study. For example, for K p = 3 to 4. 
the horizontal winds observed during the two December sol- 
stices are significantly different (solar flux levels are quite dif- 
ferent), but the two March equinox winds are similar < solar flux 
also different) and similar to the December 1991 solstice. This 
likely indicates that binning strategies for reporting mean 
w'inds need to be quite fine. Comparisons of the current data set 
with previously published data sets confirm this point. Aruhah 
et al. [1991a] binned neutral winds according to K p index and 
solar flux level, but ignored seasonal differences Conse- 
quently, even though the general trend in wind evolution is 
common to both data sets, fine details are not common, for ex- 
ample, for K p - 0 to 1. Kiruna winds demonstrate very low 
magnitude winds during dusk and magnetic midnight with 
highest amplitude near 0300 MLT for both solar minimum and 
maximum, while Watson Lake winds, during December MW 2 
solstice, are strongest at dusk and weakest post midnight and 
during March 1992 equinox weakest at dusk and strongest post 
midnight. Seasonal differences are described by Aruhah rt 
[I99lbl but unfortunately at the expense of K p binning 

The neutral wind data set described by Sica ei al | |uv M | 
was acquired during a mix of seasons and solar activity lev el- 
and is coarsely divided into three different geophysical .at iv 
ily cases. In general, the trends observed al College are mmi:!.:: 
to those observed at Watson Lake, that is. an increasing v 
stronger wind as the disturbance level increases representing a 
polar cap flow pattern. However, with the coarse College bin- 
ning it is impossible to gauge the transition period when tin* 
neutral winds convert from a midlatitudinal style How to a po- 
lar cap flow. Finally, the Watson Lake data set illustrates dial 
it is difficult to characterize neutral winds at a transitional site 
with a harmonical relationship employing diurnal, semidiur- 
nal, and higher-order terms, as was alluded to by haulun a ai 
[ I993| The combination of middle-latitude solar heating .n;d 
high- latitude momentum transfer effects on the neutral wind 
(low should certainly be diurnal I y reproducible but max n»»i 
easily be represented by tidal behavior except for the higher 
disturbance levels. In addition, there is a strong asymmetrv hr 
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tween dusk winds and dawn winds at higher disturbance levels 
which results in a better developed vortex at dusk where the ro- 
tational Coriolis force assists the ion-drag forcing rather than 
opposing it at dawn [Fuller -Rowell et ai , 1984). 

Regarding the Fabry-Pdrot interferometer scans acquired 
while serendipitously observing aurora, the potential exists to 
use the wind and temperature information to characterize mid- 
dle thermospheric altitudes. The altitude of peak emission of 
the O 1 (6300 A) airglow feature is generally believed to origi- 
nate from an altitude near or above 240 km, though the full 
width at half maximum of the emission is about 50 km [Cogger 
et al . , 1980], Auroral electron precipitation of mean energies 
in the range 1 to 10 kcV are preferentially stopped at middle 
thermospheric levels and in the process increase the plasma 
density which consequently supplies the dissociative recom- 
bination reaction producing O I (6300 A) emission. The result, 
as may be simulated by the GLOW routine, is a decrease in the 
altitude of the peak volume emission rate for this feature. It is 
characteristic for the lower altitudes to be cooler than the upper 
thermosphere and if the neutral temperature associated with the 
Fabry-Pdrot measurement is accurate, then it is possible to 
probe lower altitudes. Unfortunately, the Fabry-Pdrot on its 
own, cannot distinguish the mean energy of the precipitation 
and, more seriously, may not easily discriminate against po- 
tential contaminating molecular auroral emissions which may 
also be recorded. The contaminant emission will act to either 
enhance the recorded background level between observed 
fringes, or actually mix with the airglow channel. Further 
spectral simulation studies are required to address these prob- 
lems if true auroral observations are desired. 

Conclusions 

Observations of the thermospheric neutral winds at a transi- 
tional site near the auroral oval boundary were conducted for 
four different seasons between November 1991 and April 
1993. The neutral wind data sets were ordered based upon the 
terrestrial geomagnetic disturbance level on an hourly basis. 

1. For the lowest level of geomagnetic disturbance and the 
lowest level of solar activity, the horizontal neutral wind flow 
pattern is primarily controlled via solar heating in the torm of 
a pressure gradient force. This is most evident in the December 
1992 solstice data set [F j 0 7 = 125) tor K p = 0 to I. 

2. For the higher levels of geomagnetic disturbance, the 
neutral wind tlow pattern represents a polar cap Row driven 
primarily by momentum transfer from a rapidly moving ion 
flow. Dusk/dawn asymmetries arise from associated rotational 
forcing (Coriolis acceleration) resulting in a neutral How 
which is highly tangential and anticlockwise in the dusk sec- 
tor and less so in the dawn sector. 

3. The VSH model formulation of the NCAR TGCM repre- 
sents much of the horizontal neutral wind data correctly. Un- 
fortunately, with the coarse 5°\5° latitude/longitude grid size, 
the finer details regarding the transition between midlatitude 
flow and polar cap tlow are blended together. The correspon- 
dence between the VSH representation and the Watson Lake 
observations is best during cases of higher magnetic distur- 
bance levels when the site resides within the auroral oval. 

4. The O 1 (63(H) A) airglow signal level, in addition to be- 
ing higher when the solar flux level is greater, also increases 
as the geomagnetic disturbance level increases. I his is also 
simulated by the GLOW model Cases of distinct auroral emis- 
sion may eventually provide a diagnostic tool tor probing the 
middle thermospheric levels. 


The Fabry-P6rot interferometer experiment at Watson Lake 
has shown the potential for remotely observing the O I (6300 
A) airglow in an unattended, long-term fashion with current 
state-of-the-art hardware. A network of environmentally stable 
aeronomical observatories can provide a diagnostic tool for 
monitoring upper atmosphere thermodynamics in the same 
fashion as remote meteorological weather stations. The de- 
ployment of such a network would provide a mechanism to ob- 
serve the solar and magnetosphenc impact of "space weather" 
on the near-Earth environment and assist in the fine tuning of 
circulatory models which would aid in forecasting the same. 
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